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Abstract 

We derive expressions for the phase-space of a particle of momentum p decaying into N 
particles, that are valid for any number of dimensions. These are the imaginary parts of so- 
called 'sunset' diagrams, which we also obtain. The results are given as a series of hypergcomctric 
functions, which terminate for odd dimensions and are also well-suited for deriving the threshold 
behaviour. 

I. Introduction 



With so much attention focussed on the properties of branes embedded in higher dimensions, it 
is of interest to examine the way in which the phase space changes as the space-time is enlarged, 
since this is what primarily determines the statistical properties of multiparticle systems; only when 
CN . amplitudes are accentuated or suppressed along certain directions in space is there a pronounced 

effect on the phase space. This paper is devoted to the topic of multidimensional phase space. We 
assume that extra dimensions are spatial — though one may envisage that time-like coordinates can 
be handled by Euclidean continuation — and that all particles propagate into the bulk. Not only 
that, we suppose that the space-time is flat so that we can neglect topological aspects. If conditions 
arise that constrain the particles to some subspace, then one may obtain the reduced phase space 
by altering the number D — 1 of spatial dimensions rather trivially; or else if topological effects 
. become important we may be able to sum over discrete modes along flat directions. Thus in almost 

every respect we are simply investigating the dimensional continuation of standard 4-D expressions 
for the decay of a system carrying momentum p into any number N of particles. 

There has been much early work in this area but it is mostly in the context of four dimensions 
(or lower) and has often been confined to few particles (like two, three or four). The most pertinent 
recent results in this connection are those of Groote et al Q and of Davydychev and Smirnov || 
because they are the most general. Our own results apply to any value of N and D and are also 

singularly well-suited for deriving threshold expansions as \fp 1 ~ — > [m\ + ui2 H + tun)- Because 

phase space p is nothing but the imaginary part of the 'sunset diagram' for p — > p\ + p2 + • • -Pn, 
our procedure will be to start with this sort of diagram and then obtain pjv from its discontinuity 
in p 2 as we cross the threshold. (Pseudothreshold singularities of these diagrams exist too but are 
not relevant for physical phase space.) 

In Section II, we give the most convenient Feynman parametric description of the general 
diagram for any dimension D = 21 and show how one may readily derive the leading threshold 
behaviour by expanding about the minimum of the combined denominator; although one may in 
principle derive the next to leading behaviour by this method, and so on, we do not pursue this 
parametric approach any further, as there is a much better way of handling the problem, which is 
described in the following sections. First, in section III, we derive some special cases of the results 
(small N) by making use of phase space recurrence methods and next show how they may more 
easily be found by Fourier transformation of the propagator products in section IV. This is our cue 
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for tackling the most general situation, but because the form of our expressions makes it tricky to 
take the massless limit of any individual particle we specifically suppose that M of the particles 
are massless and N — M are massive. Our results, stated in section V, are given as a series of 
hypergeometric terms having the form (herafter F stands for the usual 2-F1 function) 

PN oc dj(.P 2 ~ <? 2 ) j F(a,b;c;l - p 2 /a 2 ), 

j=M(3/2-t)+2N(£-l)-£+l/2 

where a is the sum of the masses. An agreeable property of this expansion is that it terminates for 
odd D and is tailored for deriving the threshold behaviour. 



II. Feynman parametric form 

In 2£ dimensions, sunset diagrams with N internal lines, produce integrals of the type[], 



in(f,{»}) =i (n/ **(p-i». 

We have included a Gamma function in each numerator since the propagator powers usually arise 
via momentum or mass derivatives of the case v = 1, and it simplifies the subsequent expressions. 
By attaching a Feynman parameter a to each internal line, one may combine denominators in the 
standard fashion and integrate over the internal momenta to establish)^ that 

(_iyv (N x da , \r(£ i u i + t-M)8(l-'£ i a i ) 



= n/ 



but it is easier to prove the result by induction in fact. Suppose that (2) is true for N; the case 
N + 1 represents a further convolution: 

r N +l(p;{v})=i J ^ e kI N (k;{u})T(iy N+l )/[(k+p) 2 -m 2 N+1 }^ 

so introduce an extra Feynman parameter (3. Combining the new denominator with (2) and inte- 
grating over intermediate loop momentum k, one remains with 

r ( \ (-V N+1 (u f 1 da * \ 1\r ref^^-Wj^l-^aOll-/?)^- 1 ^- 1 
I N+1 (p) = —— I I / — — j— \ dp- 



[p 2 (l-P)-E l H/^)-m 2 N+1 (l-P)/P}^ 



Ui—Ni+v N+1 



Now all one need do is rescale a« = Pi/(1 — P) for i = 1 to N and call P = Pn+i- A final relabelling 
of P — > a reproduces (2) for N —> N + 1. Since we know the result (2) is correct for N = 2 — it is 
very familiar to graphologists as the simplest self-energy calculation — we have thereby proved the 
result inductively for any N. 

The singularities of In(p) m P 2 wm arise whenp 2 equals the combination M 2 (a) = ^^ 1 (m?/ai), 
so let us examine the behaviour of M 2 (a) in the region of integration < on < 1, subject to the 
condition J2i a i = 1- F° r definiteness, we shall assume for the rest of this section that all masses are 
nonzero. By introducing a Lagrangian multiplier for the last constraint, it is easy to see that M 2 (a) 
is minimised in the physical region when the Feynman parameters equal a^o = Tn>i/ J2jL\ m j = m i/&i 

4 The bar notation means: division by 2tv when integrating over each momentum component and a factor of 2ir 
for each delta function. 

5 This involves some finicky rescaling of the Feynman parameters. See reference M for the case N = 3. 
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whereupon M 2 (ao) = (J2jLi mj) 2 = cr 2 . Other values of CKj are possible (by reversing one or other 
of the signs of c^o) but they correspond to pseudothresholds and lie outside the integration region. 
It is then clear that the magnitude of In(p) will be dominated by «j values in the vicinity of a^, 
so it is sensible to expand about these minima if we want to determine the leading behaviour near 
threshold. 

This procedure can be developed systematically (see the Appendix for the example N = 2 with 
v \ = v i = 1) but as there is an alternative and preferable way of tackling the problem, which 
we defer to the following sections, here we shall simply extract the leading threshold behaviour, 
because this can be done rather quickly and painlessly. Begin with (2) and note that 

M 2 (a) = a 2 + ^(ccj - a i0 ) 2 cr 3 /mi - ^(aj - aiofa 4 /m 2 H 

i i 

One sees that for 1 —p 2 /a 2 = A 2 small the integral is dominated by values of a near the minimum 
«o- Therefore write a, L = c^o + A£j, so that the denominator of (2) reads 

- A V [1 + a £ g/rm + Aa 2 £ £/m? + 0(A 2 )]. 

i i 

Since the integrals in a can be made to run between and oo, because of the delta-function 
constraint, the integral over £j runs from —rrii/aA to oo. So, in the limit as A — > 0, we can take 
£j to run from — oo to oo (the correction to the integral is exponentially small as exp(— 1/A 2 )) and 
expand the product of the aj about the minimum. Assuming all masses are non-zero (see later 
sections for a relaxation of this condition) we are thereby led to the leading behaviour, 

I N (p,{u}) = c M A( JV - 1 )^ +1 )- 2 E i ^ ( J] m f-^-i)/( 47 r)^ iV - 1 ) ( 7( 2 - Ar ^- Ar+ S i ^, ( 3 ) 

i 

where 

»»-(ny*j(-D ^ [1 + (TE , 5|/mi] E,,-(->r 

If one now represents the delta function and the denominator by integrals, the coefficient c^i can 
be explicitly evaluated as follows: 



x(-l) 7V -E,^ []exp 



'(N-l)i e -a-k 2 /4a 



— 6 - tt^ 2 
nii 2aa 

= (_l) Ar ~E l "i ro V2 7r (AT-l)/2 (7 -iV/2 r( ^ v ._y + 1/2) ((A T _ w = Y[ m, 

j i 

So, all told, the leading threshold behaviour is dominated by 
iNVP,Wi) K i; (4vr)(^- 1 K a (2-N)t-N/2+^ iVi 

(4) 

The case of greatest interest is = 1, all i, when (4) reduces to 

In{p) = ^(N~l)t ^ w 3/2-/ a N/2+l(2-N) ~0- ~ P* / <**)^ N ^ {N+1)/2 . (5) 
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We may then continue this expression above threshold (p 2 > a 2 ), in order to obtain the behaviour 
of the TV-body phase space as the discontinuity: 

j^N+iyz ^-3/2 ( J pj g 2_ 1 \t{N-l)-(N+l)/2 

p N ( P ) = 2£*I N (p) ~ {47r){N _ iy r((jV -!)(£_ i/ 2) ) ( 6 ) 



and reduces to the 
when we set £ = 2, 



This result agrees with the answer obtained by Davydychev and Smirnov 
well-known four-dimensional behaviour found by Hagedorn and Almgren [jl 
namely 

PN(P) ~ Q^/ 2 , 

where Q = ^/p 2 — a is the energy release. However (6) supplies the answer in any dimension with 
the appropriate coefficient. 



III. Exact results 

We now consider phase space in general and derive some precise results for any £, N as they arise 
from recurrence relations between phase space expressions for smaller N and not necessarily around 
threshold. Begin with the definition of the iV-body phase space integral, 

p N (p) = Pp ^ 1+2+ ...N =Il(J ^ *(Pi)*(P? " m i)) ^(P " ( ? ) 

The measures d 2l ~ 1 p = \p\ 2e - 2 d\p\. (sin Of^dO^- 1 /T(l-l) = \p\ 2i - 2 d\p\.2^ 1 - 1 / 2 /T(l-l/2), come 
in useful if one were able to integrate over angles. Thus the two-body result is readily evaluated in 
this way to be 

^(4vr) 1 /2^ g 2^3 



r(£-i/2)vy ' 



where q is the centre of mass spatial momentum, so that ^ q 2 + m 2 + \J q 2 + m?, = \fp 2 ■ Tidying 
up, the result can be expressed covariantly as 



TT 1 "^ - 1)A 2 ^ 
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Pp^i+2 = 2 2^-i(p2 V )£-i r (2l- 2 ) ; A = VP 4+m i+ m 2- 2m\m\ - 2p 2 m\ - 2p 1 m\. (8) 

The three-body phase space can also be evaluated by brute force methods and reduced to a 
triangular integral over three Mandelstam variables: 

Pp^i+2+3 = o \ / / dsdtdu5(s + t + u-m\-ml-ml - p 2 )[®{s, t, u)] e ~ 2 6(®), 



{4jr) 2i r(2£ - 2) 

where $>(s,t,u) is the Kibble || cubic (simply a Gram determinant), 

j. \ j. / 2 2 , 2 2\ 2 2 , 2 2\ / 2 2 , 2 2\ 

<P(s, t, u) = stu — s(m 2 m 3 + p m 1 J — t(m 3 m 1 + p m 2 ) — u(m l m 2 + p m 3 ) 

■ n / 222, 222, 222, 2 2 2\ 2\2/|-»|2 I - |2 | -> |2\ 

J ^l\m x m 2 m z + p m 1 m 2 +p m 2 m 3 + p m^rn^ = —p A , \p2\ , |P3| J. 

By changing variables to s and t — u, one integration may be performed and the problem reduced 
to the single integral ||], 

(or, \2-U r (^2- m3 f 
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where T> = [s — (m^ + Vf 2 ) 2 ]! 5 — ( m 3 ~~ V^P 2 ) 2 ] I s ~ ( m i + m2) 2 ][s — ( m i — m 2) 2 ]- When we are in four 
dimensions (I = 2) the integral is nothing but the area within the Dalitz plot but its evaluation for 
arbitrary masses is easier said than done, because (9) is an elliptic function in general! However in 
three dimensions the integration over the region becomes possible because it can be converted into 

. . 1 f(V^-m 3 ) 2 d s if mi + m2 + m3 S 

P3{P) ' 




J(m 1 +m 2 ) 2 \fS 

In fact the three-body phase space integration is tractable for all odd D dimensions, because one 
remains with a polynomial in half-integral powers of s which is readily handled mn. Two other 
cases are amenable to an exact treatment in terms of 'elementary' functions for any £ and A" = 3, 
namely (i) two masses are set to zero, or (ii) one mass vanishes and the two other masses are equal. 
In case (i) put m,2 = = 0, mi = m so T> —* (s — p 2 ) 2 (s — m 2 ) 2 whereupon the 3-body phase 
space reads 

while in case (ii) put 1713 = 0, mi = m-2 = m so that T> — ► s(s — p 2 ) 2 (s — Am 2 ) and 

Of course we can also proceed to the limit m — * in either case. Any other set of masses produces 
'elliptic' results Jl], || for even D or integral £. 

A more systematic way of arriving at Lorentz invariant answers, without resorting to integra- 
tions over spatial momenta in the standard approach and interpreting them covariantly, is to apply 
Almgren's method. In that method one partitions the set of outgoing particles into subsets (call 
them A, B etc.) and obtains p via mass integration convolutions; for instance the three-body result 
(9) can be construed as 

Pp^l+2+3 = J Pp^A+3 PA^l+2dm 2 A , 

where s is interpreted as the intermediate m\. More generally, one may evaluate the A^-body phase 
space through the double integral 

P P ^i+2+-+j+(j+i)-+N = J Pp^A+B PA^i+2+-j Pb^(j+i)+-+n dm 2 A dm 2 B . (10) 

Although this is quite a satisfactory numerical way of calculating p, it does not shed a great deal of 
analytical light on the nature of the problem; but it does serve as a nice check of analytical results 
obtained in a different manner, which we will now exhibit. 



IV. Coordinate space method 

The sunset diagram having all v=l is nothing but the 2^-dimensional Fourier transform of the 
product of N causal functions: 

r N 

In{jP) = —i d 2e x exp(ip.x) Y[[iA c (x\mi)], (11) 
J i=i 
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and the phase space integral is simply given by 

Pp->l+2+-+N = 23f/jv(p), 

which is nonvanishing for p 2 > (m\ + 7712 + • • • + m^) 2 . Of course (11) is easier stated than done 
except in the simplest cases (like N = 2 or 3) because the causal function is a Bessel function in 
general ||, 

iA c (x\m) = * K f _i(mr); r = \/ -x 2 + ie. 

In the massless limit, when m — > 0, the sunset diagram reduces to a "superpropagator" with integer 
index since 

iA c (x|0) = iL> c (x) = T{£ - l)/4vrV 2 ^ 2 . 

However, as noted and indeed emphasized by Berends et al || and by Groote at al Q, one can 
make considerable progress using coordinate space methods in the odd-dimensional massive cases, 
because the Bessel function reduces to an exponential times a polynomial. 

Let us begin systematically by considering the most trivial case when all intermediate masses 
vanish. Carrying out the angular integration in (11), the sunset integral reduces to 

poo rrVi! 1 \ n-^ 

I N (p) -» -(2tt)V~' / dr i*J t - X {qr) 



T{i-1) 



which can be evaluated straightforwardly ]10|| . Thus one finds that 



P 2 , (12) 



j™ o i n (p) = -(4tt)^ r(iV£ _ N) (V - w > (13) 

so the massless phase space integral reduces, in 21 dimensions, to^| 

iimp KW . 29 / kW . w -"X i)) cr-^w- («) 

Both (13) and (14) are exact too. One may readily check that the cases N = 1,2,3 produce 
the correct answers, by taking appropriate limits of earliers results. Also it is very instructive to 
check that the Almgren recurrence formulae, such as (10), are properly obeyed, not only in their 
momentum dependence, but in their multiplicative coefficients. 

Now it turns out that the integral of two Bessel functions with an exponential can also be 
handled [O] . Therefore one can improve on above and calculate analytically the sunset integral for 
one particle massive (m) and the remaining N-l particles massless. In those circumstances one is 
led to |l| 

Is(p) = jf* ^ 1)M+1 ^i(*iH.(^) ri 

Noting that in respect of the variable z the hypergeometric function F(a, b; c; z) for real a, b, c has 
a branch point at z = 1 and that the discontinuity across the cut (which runs to +oo) is 

*F{a, b; c; z) = -^M^^-^^1f(c - a, c - b; c - a - b + 1; 1 - z), (16) 

r(ojr(ojr(i — a — b + c) 



6 Here and afterwards, we use the fact that OT(a)(-fc 2 + ie) a = -7r|fc 2 | a 0(fc 2 )/r(l - a) for real a. 
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we deduce the phase space result, 



Pn(p)'- 



2 2£(JV-l)-l( m 2)JV<-JV-£+i r ( 2 (JV - 1)(£ - 1)) 



F(N(£-l), (N-l)(£-l); 2(N-l) (£-!)■ 1- 



P 



to z 
(17) 

The pair of expressions (16) and (17) are precise as well. Contrast (17) with the expression found 
by Beneke and Smirnov [jl^] for a one-loop vertex diagram with two massive and one massless 
particle. By taking the limit m — > of (16) and (17) and suitably manoeuvring the hypergeometric 
function^ one can show that they collapse into (13) and (14). Additionally we may deduce the 
particular cases N = 2, 3 previously by direct substitution; for instance, in four dimensions (1 = 2), 
one gets 



Pn(p) 



2ir(p 2 — to ) 



2\2N-3 



-F(N,N- 1;2(N- 1);1 - E-)e(p 2 -m 2 ), 



(16ir 2 ) N - 1 (m 2 ) N - 1 r(2N 



m- 



so the three-body phase space with a single massive particle (p 2 > to 2 ) brings in a logarithmic 
function: 



(p 2 



TO 



2\3 



7687T 3 TO 4 



-F(3,2;4;l 



EL 

TO 2 



TO 



256vr 3 p 2 



(4-i)- 2 4m4" 

TO TO m/ 



V. Hypergeometric Expansions 



As soon as we have at least two massive particles, we come across integrals involving products of 
three Bessel functions with different arguments and a power of r. This is not given in the standard 
texts0, though it is surely some generalization of a hypergeometric function for we definitely know 
how to evaluate the N = 2 case directly in momentum space; but for larger iV it would seem that 
the coordinate space method stalls. Only when p = 0, which corresponds to a vacuum 'watermelon 
diagram', can the integrations sometimes be carried out []D|. However, we will now present a 
method which is still suitable for handling the general problem and which works wonderfully well 
in odd-dimensional spaces. It relies on the observation that the modified Bessel function possesses 
an asymptotic representation (/i = 4i/ 2 ), 



K v {z) 



1 + 



fj, - 1 ( M -l)(/x-9) + (/x-l)(/x-9)(/i-25) + ; 



8 2 



2!(8z) 5 



3!(8z) s 



that terminates for v or £ half-integral, corresponding to D odd. Even for integral £ or even 
dimensions, we shall see that it provides a very nice threshold expansion in terms of hypergeometric 



functions [14]. To show how this works, assume first that all the intermediate particles are massive. 
Since 

the product of causal Green functions produces the series 



N 



c (x\mj; 



7T 



1/2-i 



N N 



2^+l/2 r ^-l/2 



ne-3/2 - 
m 3 



i + 



(21- l)(2£-3) 



8r 



N 1 

E — 



+ ••■ 



7 e.g. using the relation F(a, b; c; z) = (1 — «) b _F(c — a, 6; c; z/(z — 1)). 



The formula 2f~r»+' K»{ar)K»{br)Mcr)dr =^/f (^)^ +1 Ig±0tiM|^ F (i/2- / ,, l/2+ M ; 3/2+^; (1-^/2) 



where 2 = (a 2 + 6 2 + c 2 )/2afe, is given in [ H . 
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The leading (and in 3-D the only) term of the rhs above yields the leading behaviour of the sunset 
diagram, 



'■Wf Jo 

2 i-Ne-N/2 7r e-Ne+N/2 T ^2-N)£ + N/2) 2£-N£+N/2 2£-N£+N/2 + l p 2 
~ T ^ a (2-N)e+N/2 w 3/2-e F ( 2 ' 2 ■J^AW) 

where w = X\ i m,i^ a = J2i m i- The imaginary part then produces the leading contribution to the 
phase space integral, 

o 2 2l(l-N)+l 7T i{l-N)+{N+l)/2( p 2_ (T 2y(N~l)~{N+l)/2 ^ ^ 

Pn{p) = ^-e^Ne-i-N^^N _ i) _ (AT + l)/2) ° ip ~ ° ) * 

It should be noted that expressions (19) and (20) represent the complete results for 3-D, when they 
reduce respectively to 

7r 3 / 2_Jv r(3 - N) 3 - N 4 -TV 3 p 2 
n(J>) ~* ~ 2 2N ~ 2 <j*~ n ^ ' ~2~' 2' ^' 

and 

N N-l AT p 2 , ^-a) N ' 2 e(p 2 -a 2 ) 



2 3 ^"V v ~ 1 r(iV - 1) v 2' 2 ' ' a 27 (47r) JV - 1 r(iV - lJVF ' 

Nonleading behaviours (for values of £ ^ 3/2) of /at and p^r can be found from expansion (18); 
for instance to the next order we encounter the terms, 

! (21 - l)(2l - 3)^« 1 - Af ) +Af / 2 r(2^ - l - JVl + jV/2) 

~~ 22+W+Af/2 ro 3/2-£ (T £(2-Ar)+Ar/2-lp^^ X 

F( U-™ + W-\ ^-^ + N/2 ](] ^ (21) 

(g£ _ 1)( 2 £ - 3) ^(l-jV)+(^V+l)/2 (;p 2 _ a 2y(N-l) + (l-N)/2 ^ 

Pn(P) - 2 2 ^-i)+3 ro 3/ 2 -V^-V2) M r(£(Af-l) + (3-iV)/2) {P a)X 

F( m ± ^N_ 1 -N^, AN _ 1) + ^ l _^ (22) 

where 1/fi = J2j ^/ m j', an d so on for higher I^,p k N . The method is entirely systematic and one 
may proceed to as high an order k as needed. Observe that the resulting series is a sort of threshold 
expansion because the powers of (p 2 — a 2 ) which multiply F(a, b; c; 1 — p 2 /o~ 2 ) keep on increasing 
as we raise k while the denominators are associated with sums of the type 2~2j(f/rrij) k ■ 

However these expansions are deficient in one respect: it is very tricky to consider the limit 
as one or several masses vanishes, because the approximation (18) is fairly useless for massless 
propagators; in that case we should be using directly iD c (x) oc r 2 ~ 2i , rather than the asymptotic 
expansion (18). Finally then we will consider the case where M of the N particles are massive 
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while the remaining N — M are massless; this covers essentially all cases of interest. The sunset 
integral is given ab initio by 



( 27 £ ir , ,fr(£-i)\ N ~ M f(7Tr) l / 2 - e \ M 



In(p) = "^r I" drr^J^qr) 
q Jo 



(2e-l)(2£-3) 



(23) 



where the symbols now refer simply to the massive intermediate particles; therefore w = YlfLi m j, 
a = J2jL\mji = X/^ii(l/ m j)- The integrals in (23) are readily performed and the leading 
term of the full answer is 

tQ ^ 2 2(M-N)+l-Ml~M/2 7T e(l~N)+M/2^^_ 1 ^N-M r ^ 2 _ M ) + ( M _ N ^ 2 £ _ 2 ) + M/2) 

zu 3/2-e a e(2-M)+(M-N)(2£-2)+M/2Y(£j 

<( 2-M) + M/2 + (M _ m ? 

Z Z 



r N (p) 



Taking its discontinuity, the leading phase space behaviour is (c = M(3/2 — t) + IN (I — l) — £ + l/2 
and p 2 > a 2 below), 

2 2 ^ 1 - Ar )+ 1 7T^ 1 - 7V ) + ( M+1 )/ 2 [T(£ - l)] N ~ M (p 2 - a 2)M(3/2-£)+2N(£-l)-(£+l/2) 
Pn(p) = zu Z/2-l (J M{7/2-Zt)+2{M+N)(£-l)-lY^ X 

F(( ,_ 1)(W _M) + M^, (< _ 1)(w _ M)+ ^i)_l ;c;1 _g ) . (25) 

All the previous cases fall out of (24) and (25) by making the relevant substitutions for M and £. 
Of course one may also derive nonleading terms 1^, p% in exactly the same way as before and they 
correspond to higher order threshold corrections. Again, these corrections entrain higher powers of 
(p 2 — a 2 ) and terminate for half-integral i. 



VI. Conclusions 

Before we can comprehend the statistical effects of phase space on brane physics, it is vital to 
understand phase space in fiat D-dimensional spacetime for any number of particles, with arbitrary 
masses. This paper has been devoted to that subject and we have arrived at results for pn(p) and 
In{p) that have culminated in formulae (14), (23), (24) and (25). These comprise the high-energy 
and low-energy characteristics and at any energy in-between. We believe that these analytical 
expressions are as compact as one can make them and will turn out to be practically useful. 
Otherwise one will be obliged to resort to numerical methods. 
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Appendix. At the threshold of a sunset 

In this appendix we show explicitly how the threshold expansion can be carried out in momentum 
space via Feynman parameters. We illustrate the case N = 2 for any dimension 21 to keep the 
algebra simple, when the sunset integral is just 



{-AttY Jo Jo \p 2 aia 2 - m\a 2 - rn\ai} 2 r 

Never mind the fact that I 2 can be expressed as a linear combination of two hypergeometric 
functions (a = I — 1, b = 2 — £, c = £) with arguments |[(mf — m 2 ±p 2 )/\ + 1] or that it can 
be accorded an exact geometrical interpretation our purpose here is to show how one may 
arrive at a systematic threshold expansion, by expanding about the minimum of the denominator, 
regarded as a function of a. Thus change variables to 

oti = — + A Pl , p 2 = (m 1 + m 2 ) 2 [l- A 2 ]. 

mi + 1^2 

The delta-function constraint on a corresponds to p\ = —p 2 = p, whereupon the sunset integral 
reduces to 

r(2 — £) .2i-3 f (™i+™2)^ dp 



hip) = ^r^A 



J- (mi 7l 2) A [™i m 2 + P 2 {mi + m 2 ) 2 (l - A 2 ) - Ap(m 2 - m 2 ,)] 

The leading behaviour is found by setting A = in the denominator of the integrand and extending 
the limits to ±oo. Since 

dp _ (m 1 m 2 ) £ " 3/2 f°° dw" 1 / 2 _ Kms)"/ 2 y/jrT(3/2 - £) 



lo [mim 2 +p 2 (mi + m 2 ) 2 ] 2 1 2{m\ + m 2 ) Jo {l+v) 2 1 2{m\ + m 2 ) T{2 
it follows that 



(m^A 2 )^ 2 ^ 
I2{P) ~ (4^(m 1 + m 2 ) r(3/2 l) ' 



However we are keen to obtain all the nonleading terms and this can be accomplished by 
rewriting the sunset integral as 

r(2-Q( mi m 2 A 2 )"/ 2 f- +j 2 e _ 2 

h(p) = . ^ / du[l + u -urj\ , 



where the new integration variable is u = py / p 2 /mim 2 , the limits are u + = / ™* v, J p and 



mi 



1 ->- w p , and the new expansion variable is ri = A ^, m ' 1 . The series in rj reads 

Tf . ( mi m 2 A 2 Y- 3 / 2 ^ T(2 + n-£) njj TJ m e _ 2 _ n 

(4vr)VP 2 „ =0 n] Ju - 

In evaluating the coefficients U n , we make use of the fact that the limits are large and we should 
distinguish between the cases of n even and n odd. For n even, we split the integral into 

U n = |2/ - / - / ] duu n (l + u 2 f- 2 - n 

OO /"OO \ 



'|u_ J\u+\J 

r((n + l)/2)r((3 + n)/2-£) 



r(2 + n-^) W| U _| J\ u+ \ 



/■OO roo \ 

+ du/- 4 -"(l + l/« 2 ) 

U_ I J In J. I / 
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where we may expand the latter two integrals in powers of 1/u 2 — which provides a further sub- 
expansion in powers of A. Thus we get 

r(2 + n-^)C/ n = r((n + l/2)r((3 + n)/2-£) + ^ 1 ' ' ' - 



k=0 



k\{2t - 3 - n - 2k) 



| u |2£-3-n-2fc _|_ 1 |2£-3-n-2fc 

On the other hand, the case of odd n can be treated directly by expanding the integrand below in 
powers of 1/v: 



u n = \ r l dw i -^ +r ^ 2 (i+\/vf- 2 - n 



2 Jul 

00 , \k 



so 



i(2 + n i)Un ~ ^ k\(2i-3-n-2k) 



1 i2l-3-n-2k _ 1 |2£-3-rt-2fc 



Again one encounters a sub-series in A. (The two series make up the 'h-h' and 'p-p' contributions 
of ref ||] at one stroke.) In short we see that the sunset integral, and of necessity its imaginary 
part, can be systematically expanded as a series in A which is proportional to the Q-value of the 
reaction by appropriately handling the Feynman parametric representation. It is possible to treat 
the sunset diagram with iV intermediate particles in a similar way, but the method becomes rather 
unwieldy, which is why we turned to coordinate space methods in sections IV, V. 
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